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Abstract 

We study the phase transitions of N-components generahzed hehmagnets 
that obey the symmetry breaking pattern O(A^) x 0(A^ — 1) ^ 0(A^ — l)diag- 
In the neighborhood of two dimensions, a = 2 + e renormaUzation group 
study reveals a rich fixed point structure as weh as a nematic-hke phase with 
partial spin ordering. In the physical case 0(3) x 0(2) — > 0(2)diag, rele- 
vant to real magnets with noncollinear ordering, we show that this implies 
an XY-like transition between the ordered phase and the nematic-like phase. 
A non-Abelian mean-field calculation qualitatively valid above four dimen- 
sions is shown to lead to the same picture but then the principal chiral fixed 
point, which had been proposed earlier as the relevant fixed point for D=3 

helimagnets, plays no role due to the appearance of a first-order line. 
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Several physical systems display a peculiar critical behaviour associated with helical 
or noncoUinear ordering. These include helimagnets, the dipole-locked phase of superfiuid 
Helium as well as Josephson junction arrays in a transverse magnetic field |I|]. A prototypical 
example is the stacked triangular lattice antiferromagnet (STA) whose low-temperature 
ordered phase breaks completely the rotation symmetry group. When generalized to N- 
component spins, the Landau- Ginzburg free energy for this system contains two quartic 
invariants: 

+U (01 + 0'^^'^ + V ({01 ■ 02f - 01 0t) ■ (1) 

When M,f > 0, this free energy describes a second-order phase transition between a high- 
temperature phase which is O(A^) x 0(2)-symmetric and a low-temperature phase with 
lower symmetry 0(A^ — 2) x 0(2)diag- This theory has been studied in the neighborhood 
of four dimensions by the standard e-expansion The main point of these studies is 

that when the number of components is large enough there is a fixed point with u*,v* 
that describes a critical behaviour different from the well-studied behaviour of the N-vector 
model. In the neighborhood of the upper critical dimension, D=4, there is a dividing 
line N(.{D) = 21.8 — 23. 4e + O(e^) in the (N, D)-plane above which there is a second- 
order phase transition and below which there is a fluctuation-induced flrst-order transition 
(no stable flxed point). This is similar to the case of the normal-superconducting phase 
transition [Q. The fate of the line Nc{D) is not known outside the e-expansion. Numerical 
studies on the STA lattice (thus in D=3) have shown that the ordering transition is second- 
order with exponents that do not belong to the O(A^) Wilson-Fisher universality classes 
u = 0.585(9), ^jv = 2.011(14) These exponents are also observed in the body- 
centered tetragonal antiferromagnet p, suggestive of a new universality class, usually called 
chiral universality class. It has been suggested 0,^,0 that these exponents are ruled by the 
nontrivial flxed point which exists for > Nc{D). 

The study of these phenomena has been also performed from the low-temperature expan- 
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sion. A non-linear sigma model that captures the Goldstone modes of the theory has been 
constructed by Dombre and Read |T^. It involves an order parameter which is a rotation 
matrix instead of a vector as in usual magnetic systems. This is due to the fact that the 
rotation group is fully broken in the low-temperature phase of noncollinear magnets. The 
symmetry breaking pattern is S0(3) x S0(2) S0(2)diag where the internal S0(2) rotation 
acts upon the 1,2 indices of the fields of Eq.(|l|). We can parametrize the rotation matrix 
of the non-linear sigma model by three orthogonal unit vectors 0i, 02, 03 and the Euclidean 
action can be then written as: 

S = fd^x^ ((V0i)2 + (V02)') + :5^(V03)'. (2) 

J iQi ^ zg2 

The STA corresponds to a bare value l/g2 = 0. The renormalization group flow of the two 
couplings gi,g2 has been studied in a D = 2 + e expansion. It has been noted that 
a remarkable symmetry enhancement happens on the line gi = g2- In this case the global 
symmetry of Eq. (Q) is SO (3) x SO (3) broken down to the diagonal subgroup SO (3). This 
is the so-called principal chiral model which is renormalizable and thus the peculiar line 
91 = 92 is stable under the RG flow. Since S0(3) x S0(3) = S0(4), this line is described 
by the usual SO (4) sigma model, at least within the D = 2 + e expansion. Since this fixed 
point is stable and has an exponent u ~ 0.74 there is a clear conflict with numerical and 
experimental findings in D = 3. 

In this Letter we shed new light on this problem by studying a new generalization to 
N components of this sigma model as well as by a non-Abelian mean-field calculation. We 
show that the sigma model has a much richer structure than previously expected and we 
propose a scenario in which the principal chiral fixed point with SO (4) symmetry plays no 
role in D = 3 due to the appearance of a first-order transition. 



In the two-dimensional case, the perturbative beta functions |12| , p!l| indicate a flow which 
is infrared unstable away from the origin gi = g2 = but the whole line = is a 
line of unstable fixed points as can read from the RG formulas. This fact has a simple 
interpretation: when g2 ^ the vector 03 becomes frozen and the remaining degrees of 



freedom are simply the rotations of the vectors (f)i,(j)2 in the plane orthogonal to (p^- Since 
(pi and 02 are themselves orthogonal, we are left with an XY model whose coupling is gi. 
The corresponding beta function is zero in perturbation theory. However it is important to 
note that in D = 2 this XY model will undergo the Kosterlitz-Thouless (KT) transition for 
a finite value of gi. This phenomenon is not seen in perturbation theory and is not found in 
the perturbative beta functions. We expect that the line of fixed points (72 = will end at 
some critical value. Above this value, the flow will go to the high-temperature fixed point 
in the disordered phase of the XY model. This transition is the unbinding transition of the 
vortices of the XY model. If we unfreeze the vector 03 by setting (72 7^ 0, then topological 
defects survive because, as noted by Kawamura and Miyashita [0, ni(S0(3)) = Z2. All 



XY vortices with even winding number become topologically trivial while those with odd 
winding number become all equivalent and nontrivial. 

To complete the phase diagram in the [gi — (72)-plane, we note that for gi = 00 the 
model d^) becomes the 0(3) sigma model which has no transition in D=2: the coupling g2 
flows continuously to the high-temperature fixed point. If we assume that there are no other 
fixed points, we are led to propose the phase diagram shown in Fig. 1. This means that 
for nonzero values of the couplings the model is always disordered and there is no phase 
transition at nonzero temperature as expected for a system with non-abelian symmetry. 

We note that the KT point on the gi axis cannot be studied by extending the model 
simply to components. The low-temperature limit of the theory (1) is an O(A^) x 0(2) 



sigma model |TJ] that has also two coupling constants. The corresponding phase diagram 
is topologically similar to Fig. 1: the KT point remains at finite coupling unreachable by 
perturbation theory. 

To clarify the situation, we will study another generalization of the 0(3) x 0(2) model 
which involves A^ orthogonal unit vectors with A^ components (pi, ... , (p^- and the action is: 

S=f rf^X-^ ((V0i)2+ . . . + (V0^_i)^) + ;r^(V0^)l (3) 
J igi ^ ' zg2 

The symmetry breaking pattern is now SO(A^) x SO(A^ — 1) ^ SO(A^ — l)diag- This is a 



sigma model which is defined on a space which is homogeneous but not maximally symmetric. 
When Qi = g2, one is dealing with the sigma model defined by the maximally symmetric 
space SO(A^) x SO(A^)/SO(A^)diag, this is the so-called principal chiral model. When g2 — 
0, it reduces to the principal chiral SO(A^ — 1) model, and when (71 ^ oo to the O(A^) 
vector model. In fact, this model has been studied some time ago by Friedan who has 



shown that it is renormalizable in two dimensions with only two coupling constants. The 
corresponding RG fiow in = 2 + e is depicted in Fig. 2. 

The principal chiral models are known to have a fixed point within the D = 2 + e 
expansion: these are the points Cn and Cat-i in Fig. 2. The O(A^) vector model has also 
a fixed point On which is at distance e from the upper left-hand corner of the Fig. 2. 
The novelty is that we find a fixed point P n with nontrivial values of the couplings gi , g2 
which is not on the diagonal gi = g2- This fixed point has two directions of instability 
and thus there are two phase transition lines: one that goes from Pat to On and one that 
goes from Pjv to Cjv-i- This implies that there is an intermediate phase in addition to the 
high-temperature paramagnetic phase and the low-temperature ordered phase. In this new 
phase, the vector 0Ar is ordered because it is in the low-temperature regime of the O(A^) 
model but the remaining — 1 vectors 0i, . . . ,4>n-i are still fiuctuating in the subspace 
orthogonal to (pN- Due to this partial ordering, we will refer to this phase as "nematic" . The 
transition between the fully ordered phase and the nematic phase is governed by the fixed 
point Civ-i and the transition between nematic and paramagnetic phases by the O(A^) fixed 
point. The fixed point Cn which has the highest symmetry governs the critical behaviour 
of the direct transition between the paramagnetic phase and the fully ordered phase. 

This phase diagram thus points to a possible explanation of what happens in the O(A^) x 
0(2) sigma model above two dimensions. Here the point Cat-i is replaced by an 0(2) point 
which is not seen in the D = 2 + e expansion. This point is located at (72 = but gi = 0(1) 
and thus the perturbative fiow of the coupling gi is that of a low-temperature phase. The 
0(3) point is seen in perturbation theory [|r^. The mixed fixed point will also be invisible 
in perturbation theory. In fact, in the fiow equations of the generalized model, the limit 



^ 3 leads to (7* — oo and = 0{e). When D ^ 2 the fixed points that were of order 
e collapse to the origin but the 0(2) fixed point becomes the KT point at the end of a line 
of perturbative fixed points. So the picture we have constructed leads naturally to a phase 
diagram as Fig. 1 for all O(A^) x 0(2) models. 

These D = 2 + e expansion results clearly show that the phase diagrams of the these 
sigma models is much more complex than previously thought Notably the fixed point 
Pat is a prominent candidate to interact with the principal chiral point Cat. To deepen 
our understanding, we now turn to a mean- field study of the 0(3) x 0(2) sigma model 
on a hypercubic lattice. This lattice model is regularized both in the ultraviolet and in 
the infrared limit: this is a sensible way to define the sigma model beyond perturbation 
theory. Note that to construct a continuum field theory one has to find a second-order 
phase transition which is not guaranteed a priori. We introduce a field Aj conjugate to each 
vector (pi in the action (^. The standard mean- field method then involves a non-Abelian 
integration over the 0(3) matrix (0i,02;03)- We start from the Hamiltonian: 

H = - KiM^)-My) = -T.K^^^-^''-$^^ (4) 

<x,y> 

where the sum is over nearest-neighbor lattice sites, Vxy is the connectivity matrix and we 
take only two distinct couplings: Ki ~ 1/gi for 0i,02 and K2 ~ 1/(72 for 03. We write the 
partition function as: 

Z = Jd\, e-*"^'-^"'"^^ J t//i(0,) e?'''^\ (5) 

Here dfi is the Haar measure on SO (3) and T is the temperature. The mean- field theory is 
obtained by a saddle-point treatment of the integral over the auxiliary fields Aj. We have 
searched spatially uniform solutions since the model (|^) is ferromagnetic. There is a phase 
where all the expectation values (Aj) are zero: this is the high-temperature fully disordered 
phase and there is also a phase where all the (Aj) are non zero: this is fully ordered phase 
that breaks the full rotational invariance. But we also find a phase where (A3) 7^ but 
(-^1,2) = 0: this is the nematic phase. The corresponding phase diagram is sketched in 
Fig. 3. 
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The most remarkable feature is the appearance of a first-order hne that crosses the 
diagonal Ki = K2. Its existence is simple to understand: the non-Abelian integral in Eq.(^ 
has no simple closed form but its expansion in powers of group invariants built from the 
matrix A = (Ai, A2, A3) is simple. Since we are dealing with the rotation group there is 
an odd invariant that appears in the mean-field potential: Det[A]. This cubic contribution 
gives rise to the first-order transition (thick line in Fig. 3). This line terminates at tricritical 
points and is continued by second-order transition lines: in particular for K2 small enough 
there is a second-order transition to the fully ordered phase. In the limiting case K2 = 0, it 
is not necessary to introduce the auxiliary field A3 in Eq. (|^) and a standard calculation leads 
immediately to the Landau theory in Eq.(|l]). When < K2 < Ki, the field A3 is necessarily 
present but it remains massive at the phase transition: At the transition the fields Ai, A2 get 
a non-zero expectation value, but the det[A] term in the potential acts as a magnetic field 
and immediately induces also an ordering of A3. 

The intermediate ordering transitions that occur in the K2 > Ki region are second-order 
close to the boundaries. These lines still exist in the unphysical region below the first-order 
line: they are then spinodal lines and they converge right at the diagonal toward the chiral 
point C3, which is thus metastable. As a consequence, we note that the latticized principal 
0(3) chiral model does not lead to a continuous theory: there is no place where the correla- 



tion length diverges. In fact there is numerical evidence from Monte-Carlo studies |lT6| , p!7[| for 
a first-order transition in the model (|) at Ki = K2 in three dimensions in agreement with 
the mean-field prediction. These studies have also obtained marginal evidence for the chiral 
universality class exponents of Refs.( [^0]) by simulation of (^) for the value K2 = 0. This 
Hamiltonian is expected to be in the same universality class as the STA-type helimagnets. 
From our findings it is however clear that the proximity of a tricritical point as seen in 
mean-field theory may lead to difficulties in the observation of the true critical behaviour: 
indeed the first-order line begins at K2/K1 = 8.5. 

In conclusion, we have shown the existence of a nematic phase with partial spin ordering 
in the family of sigma models SO(A^) x SO(A^ - 1) ^ SO(A^ - l)diag- For > 3, all relevant 



fixed points are captured hj a, D — 2 + e expansion. In the physical case SO(3)xSO(2) 
there is an XY phase transition between the fully ordered phase and the nematic phase. We 
obtain a similar picture from mean-field theory with the appearance of a first-order line that 
continue the helimagnetic second-order line between the paramagnetic and the fully ordered 
phases, and isolates the principal chiral fixed point C3 with SO(3)xSO(3) symmetry in the 
metastability region. The simplest scenario is that this line appears at an unknown critical 
dimension Dc above which the D = 2 -|- e is replaced by the mean-field picture. In this 
respect we note that the intersection of the spinodal lines on the diagonal is suggestive of a 
collapse of the two fixed points Pn and Cn- If this critical dimension Dc is between 2 and 
3, there is a natural explanation to the fact that the chiral universality class has exponents 
different from 0(4). 
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FIGURES 
0(3) 




FIG. 1. Global renormalization group flow of the 0(3) x 0(2) sigma model in two dimensions. 
The perturbative line of fixed points on the gi axis ends at a KT fixed point. 
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FIG. 2. The RG flow and the corresponding phase diagram of the 0{N) x 0{N — 1) sigma 
model in the neighborhood of two dimensions D = 2+e. There is a high-temperature paramagnetic 
phase, a fully ordered phase at low temperature and an intermediate nematic-like phase where the 
4>N vector is ordered but the other vectors are disordered. The diagonal gi = g2 is the principal 
chiral model. 
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FIG. 3. The non-Abelian mean-field phase diagram of the 0(3) x 0(2) sigma model (sketchy). 
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